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Abstract 


ATLAS has measured two-particle correlations as a function of relative azimuthal-angle, A0, 
and pseudorapidity, A?/, in V^=13 and 2.76 TeV pp collisions at the LHC using charged par¬ 
ticles measured in the pseudorapidity interval |?7|<2.5. The correlation functions evaluated in 
different intervals of measured charged-particle multiplicity show a multiplicity-dependent 
enhancement at A0 ~ 0 that extends over a wide range of A?/, which has been referred to 
as the “ridge”. Per-trigger-particle yields, Y{A(p), are measured over 2<|A?7|<5. For both 
collision energies, the Y{A(p) distribution in all multiplicity intervals is found to be consis¬ 
tent with a linear combination of the per-trigger-particle yields measured in collisions with 
less than 20 reconstructed tracks, and a constant combinatoric contribution modulated by 
cos {2A(p). Tbe fitted Fourier coefficient, V 2 p, exhibits factorization, suggesting that the ridge 
results from per-event cos {2<p) modulation of the single-particle distribution with Fourier 
coefficients V 2 - The V 2 vaiues are presented as a function of multiplicity and transverse mo¬ 
mentum. They are found to be approximately constant as a function of multiplicity and to 
have a pj dependence similar to that measured in p-i-Pb and Pb-i-Pb collisions. The V 2 values 
in the 13 and 2.76 TeV data are consistent within uncertainties. These results suggest that the 
ridge inpp collisions arises from the same or similar underlying physics as observed in p-i-Pb 
collisions, and that the dynamics responsible for the ridge has no strong ^/s dependence. 
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Measurements of two-particle angular correlations in high-multiplicity proton-proton {pp) collisions at 
a center-of-mass energy ^[s = 1 TeV at the LHC showed an enhancement in the production of pairs 
at small azimuthal-angle separation, A 0 , that extends over a wide range of pseudorapidity differences, 
A 77 , and which is often referred to as the “ridge” [1]. The ridge has also been observed in proton-lead 
(p-tPb) collisions [2-7], where it is found to result from a global sinusoidal modulation of the per-event 
single-particle azimuthal angle distributions [3-6]. While many theoretical interpretations of the ridge, 
including those based on hydrodynamics [8-12], saturation [13-23], or other mechanisms [24-30], have 
been, or could be applied to both pp and p-i-Pb collisions, it has not yet been demonstrated that the ridge 
in pp collisions results from single-particle azimuthal anisotropies. Testing whether the ridges in pp and 
p-i-Pb collisions arise from the same underlying features of the single-particle distributions may provide 
insight into the physics responsible for the phenomena. Separately, a study of the ^[s dependence of the 
ridge in pp collisions may help distinguish between competing explanations. 

This letter uses 14 nb“^ of = 13 TeV data and 4.0 pb“' o^ ^ - 2.76 TeV data recorded during LHC 
Run 2 and Run 1, respectively, to address these issues. The maximum number of inelastic interactions per 
crossing was 0.04 and 0.5 for the 13 and 2.76 TeV data, respectively. Two-particle angular correlations are 
measured as a function of A 77 and A0 in different intervals of the measured charged-particle multiplicity 
and different px intervals spanning 0.3<pt<5 GeV: 0.3-0.5 GeV, 0.5-1 GeV, 1-2 GeV, 2-3 GeV, 3- 
5 GeV. Separate px-integrated results use 0.5<px<5 GeV. Per-trigger-particle yields are obtained from 
the long-range (|A? 7 |> 2 ) component of the correlation. A new template-fitting method is applied to these 
yields to test for sinusoidal modulation similar to that observed in p-i-Pb collisions. 

The measurements were performed using the ATLAS inner detector (ID), minimum-bias trigger scintil¬ 
lators (MBTS), forward calorimeter (FCal), and the trigger and data acquisition systems [31]. The ID 
detects charged particles within |? 7 l< 2.5 using a combination of silicon pixel detectors, silicon micro-strip 
detectors (SCT), and a straw-tube transition radiation tracker (TRT), all immersed in a 2 T axial magnetic 
field [32, 33]. The MBTS sysfem defecfs charged parficles using fwo hodoscopes of counfers positioned 
ntz- ±3.6 m. The FCal covers 3 . 1 <|? 7 |< 4.9 and uses fungsfen and copper absorbers wifh liquid argon as 
fhe acfive medium. Befween Run 1 and Run 2, an additional, innermosf pixel layer was added fo fhe ID 
and fhe MBTS was replaced. 

The ATLAS trigger system [34] consists of a Level- 1 (LI) trigger implemented using a combination of 
dedicated electronics and programmable logic, and a software-based high-level trigger (HLT). Charged- 
particle tracks were reconstructed in the HLT using methods similar to those applied in the offline analysis, 
allowing triggers that select on the number of tracks with px>0.4 GeV associated with a single vertex. For 
the 13 TeV measurements, a minimum-bias LI trigger required one or more signals in the MBTS while 
the high-multiplicity trigger (HMT) required at least 900 SCT hits and at least 60 HLT-reconstructed 
tracks. For the 2.76 TeV data the minimum-bias trigger selected random crossings at LI and applied a 
threshold to the number of SCT and pixel hits in the HLT, while several HMT triggers were formed by 
applying thresholds on the total FCal transverse energy at LI and different thresholds on the number of 
HLT-reconstructed tracks. HMT triggers are only used where their multiplicity selection is more than 
90% efficient. The inefficiency of the HMT triggers does not affect the measurements presented in this 
paper. This has been checked by comparing the results obtained with and without the HMT-triggered 
events, over the range where the HMT is not fully efficient. 

Charged-particle tracks and collision vertices are reconstructed in the ID using algorithms that were re¬ 
optimized between LHC Runs 1 and 2 [35]. Tracks used in the analysis are required to have px>0.3 GeV, 
| 7 /|< 2.5 and to satisfy additional selection criteria that differ slighfly between the 2.76 [4] and 13 TeV [36] 
data. 
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Events used in the analysis are required to have at least one reconstructed vertex. For events containing 
multiple vertices (pileup), only tracks associated with the vertex having the largest ^ where the sum 
is over all tracks associated with the vertex, are used. The measured charged-particle multiplicity, is 
defined as the number of tracks having pt> 0.4 GeV associated with this vertex. The distributions of 
are shown in Fig. 1. The structures in the distributions result from the different HMT trigger thresholds. 



Figure 1: Distributions of the multiplicity, of reconstructed charged particles having pt> 0.4 GeV for the 2.76 
(left) and 13 TeV (right) data used in this analysis. 


The efficiency, e{pT, rf), of the track reconstruction and track selection requirements is evaluated us¬ 
ing simulated non-diffractive pp events obtained from the pythiaS [37] event generator (A2 tune [38], 
MSTW2008FO PDFs [39]) that are passed through a GEANT4 [40] simulation of the ATFAS detector 
response and reconstructed using the algorithms applied to the data [41]. The efficiencies for the two 
data sets are similar, but differ due to changes in the detector and reconstruction algorithms between 
Runs 1 and 2. In the simulated events, the efficiency reduces fhe measured multiplicify relative to the 
pythia8 pt>0.4 GeV charged-particle multiplicity by approximately multiplicity-independent factors of 
1.18 ± 0.05 and 1.22 ± 0.05 for 13 and 2.76 TeV data, respectively. The uncertainties in these factors 
result from systematic uncertainties in the tracking efficiencies, which are described in detail in Ref. [36]. 
Those systematic uncertainties vary with pseudorapidity between 1.1% (central) and 6.5% (forward) and 
result from uncertainties on the material description. 

The present analysis follows methods used in previous ATFAS two-particle correlation measurements 
in Pb-i-Pb and p-i-Pb collisions [4, 6, 42-44]. Two-particle correlations for charged particle pairs with 
transverse momenta p^ and are measured as a function of A0 = - 0*’ and = rf - with 

|A77|<5, determined by the acceptance of the ID. The particles a and b are conventionally referred to as 
the “trigger” and “associated” particles, respectively. The correlation function is defined as: 


C{Ap,A(p) 


S{Ap,Acf>) 
B{Ap, A(p) ’ 


( 1 ) 


where S and B represent the same event and “mixed event” pair distributions respectively [45] . When 
constructing S and B, pairs are weighted by the inverse product of their reconstruction efficiencies 
l/(e(p^, p^)e{p'^, p'°)). Detector acceptance effects largely cancel in the S/R ratio. 
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Figure 2; Two-particle correlation functions, C(Ari,A^), in 13 TeV pp collisions in intervals 0-20 (left) and 
> 120 (right) for charged particles having 0.5<p^'*’<5 GeV. The distributions have been truncated to suppress the 
peak at Arj-A(p-Q and are shown over |?7l<4.6 to avoid statistical fluctuations at larger IA77I. 


Examples of correlation functions in the 13 TeV data are shown in Fig. 2 for intervals 0-20 (left) 
and >120 (right), respectively, for 0.5<p^’*’<5 GeV. The C{Ar},A(p) distributions have been truncated at 
different maximum values to suppress a strong peak at Arj - Acf) - 0 that arises primarily from jets. The 
correlation functions also show a Aty-dependent enhancement centered at Atp - n, which is understood 
to result primarily from dijets. In the higher interval, a ridge is observed as the enhancement near 
A(p -0 that extends over the full At] range of the measurement. 

One-dimensional correlation functions, C{A(p), are obtained by integrating the numerator and denomi¬ 
nator of Eq. 1 over the long-range part of the correlation function, 2<|A77|<5. These are converted into 
“per-trigger-particle yields,” E(A0), according to [4, 6, 45]: 


Y{Acf>) - 


' J B{A(f>)dAf 
f dAcp ^ 


C(A0), 


( 2 ) 


where denotes the efficiency-corrected total number of trigger particles. Results are shown in Fig. 3 
for selected intervals in the 13 and 2.76 TeV data, for the ranges 0.5<p^’*’<5 GeV. Panel (a) in the 
figure shows Y{A(p) for 0< A^^“<20 for bofh collision energies; fhese exhibif a minimum af A^ = 0 and a 
broad peak af A0 ~ tt fhaf is undersfood fo resulf primarily from dijefs buf may also include confribufions 
from low-pT resonance decays and global momenfum conservation. The higher Y{A<p) values for fhe 
2.76 TeV dafa are due fo fhe relafive inefficiency of fhe 2.76 TeV friggers for fhe lowesf mulfiplicify 
evenfs, which resulfs in larger for fhe 2.76 TeV dafa in fhis inferval. Panels (b), (d) and (f) 
show resulfs from 13 TeV dafa for fhe 40-50, 60-70, and > 90 infervals, respecfively. Panels (c) and 
(e) show fhe resulfs from 2.76 TeV dafa for 50-60 and 70-80 infervals, respecfively. Wifh increasing 
fhe minimum af A0 = 0 fills in, and a peak appears and increases in amplifude. 

To separafe fhe ridge from angular correlafions presenf in low-mulfiplicify pp collisions, a femplafe tiffing 
procedure is applied fo fhe Y{A(p) disfribufions. Mofivafed by fhe peripheral subfracfion mefhod applied 
in p-i-Pb collisions [4], fhe measured Y{A<p) disfribufions are assumed fo resulf from a superposition 
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Figure 3; Per-trigger-particle yields, F(A(^), for 0.5</7“’'’<5 GeV in different intervals in 2.76 and 13 TeV data. 
Panel (a): 0<A^‘“<20 for both data sets. Panels (c) and (e): 50-60 and 70-80 intervals for 2.76 TeV data. 
Panels (b), (d) and (f): 40-50, 60-70, and >90 intervals for 13 TeV data. In panels (b)-(f), the open points and 
curves show different components of the template (see legend) that are shifted, where necessary, for presentation. 


of a “peripheral” Y{A(p) distribution, scaled up by a multiplicative factor and a constant modulated by 
cos(2A0). The resulting template fit function, 

ytempl(A^) ^ p FPeriph(A<^) + , (3) 

where 

- G (1 + 2v2,2 cos {2A(I>)) , (4) 

has two free parameters, F and V 2 ^ 2 - The coefficient, G, which represents the magnitude of the combi¬ 
natoric component of is fixed by requiring that dA^ ytempi _ y peripheral 
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distribution is obtained from the 0<A^‘^‘^<20 interval. In the fitting procedure, the is calculated ac¬ 
counting for statistical uncertainties in both Y{A<p) and distributions. 

Some results of the template fitting procedure are shown in panels (b)-(f) of Fig. 3; a complete set of fit 
results is provided in Ref. [46]. The scaled TP®’'‘p^(A 0) distributions shifted up by G are shown with open 
points; the functions shifted up by F’TP“‘p*’( 0) are shown with the dashed lines; and the full fit 

function is shown by the solid curves. The function in Eq. 3 successfully describes the measured Y(A(p) 
distributions in all intervals. In particular, it simultaneously describes the ridge, which arises from 
an interplay of the concave TP®"p*’(A(;&) and the cosine function, the height of the peak in the Y{A(p) at 
A(p ~ n, and the narrowing of that peak which results from a negative contribution of the 2v2,2 cos {2A(p) 
term in the region near A(f> = 7r/2. The agreement between the template functions and the data allows for 
no significant -dependent variation in the width of the dijet peak at A0 = tt except for that accounted 
for by the sinusoidal component of the fit function. Including additional cos (3A0) and cos (4A0) terms 
in Eq. 4 produces changes in the extracted V 2,2 values that are negligible compared to their statistical 
uncertainties. 

Previous analyses of two-particle angular correlations in pp, p-i-Pb, and Pb-i-Pb collisions have tradition¬ 
ally relied on the “zero yield at minimum” (ZYAM) hypothesis to separate the ridge from the dijet peak 
at A(p ~ 7T. In the ZYAM method, the ridge is functionally defined fo be T(A0) - Fmin over fhe resfricfed 
range |A0| < 0min> where is the location of the minimum of Y{A(p) and Fmin = F(0niin)- However, the 
F(A0) distributions measured in low-A^^‘^ bins are concave in the region near A(p ~ 0. As a result, if the 
ridge and dijet correlations add - an assumption that is implicit in all previous analyses using the ZYAM 
method and is explicit in the template method used here - then the ZYAM method will both underestimate 
the ridge yield and produce values that vary, unphysically, with the ridge amplitude. In contrast, the 
template method used here explicitly accounts for the concave shape of the peripheral Y{A(p). Thus, the 
template fitting procedure, for example, extracts a non-zero ridge amplitude from the = 2.76 TeV, 
50 < < 60 F(A0) distribution (middle left panel of Fig. 3) which is approximately flat near A(p ~ 0, 

and would, as a result, have approximately zero ridge signal using the ZYAM method. 

Previous p-i-Pb analyses used the peripheral-subtraction method, but applied the ZYAM procedure to the 
peripheral reference and, so, subtracted F(0) from FP®’’p*’(A^). Such a subtraction will necessarily change 
the V 2^2 values, and, when applied to the 13 TeV data, it reduces the measured V 2 p by a multiplicative factor 
that varies from 0.4 to 0.8 over 30<V‘^‘^<130 [46]. However, if, as suggested by the data, FP®‘^‘P*’(A(;f>) 
contains not only a hard component, F^”‘^(A0), but also a modulated soft component, 

FPeriphf A0) = + Gq (1 + 2i;° 2 COS (2A0)), (5) 

the peripheral ZYAM method will subtract 2FGqv^ 2 cos (2A0) as part of the template fit, thereby reducing 
the extracted ^ 2 , 2 - In contrast, the procedure used in this analysis subtracts EGo(l - 1 - 2^22 cos (2A0)), 
which reduces G in Eq. 4 but has less impact on ^ 2 , 2 - In particular, if ^22 equal to the real ^ 2,2 in a 
given interval, there will be no bias. Since the measured t; 2,2 is approximately -independent, the 
bias resulting from the presence of V 2,2 in the peripheral sample is expected to be small. Thus, the use 
of the non-subtracted peripheral reference is preferred over the more strongly biased ZYAM-subtracted 
reference. 

If the cos(2A0) dependence of Y{A<p) arises from modulation of the single-particle cp distributions, then 
V 2,2 should factorize such that V 2 p{.p^,p\) = V 2 {p^)v 2 {p\) [42-44], where V 2 is the cos(20) Fourier coef¬ 
ficient of the single-particle anisotropy. To test this, the analysis was performed using three intervals: 
0.5-5, 0.5-1, and 2-3 GeV with 0.5<p^<5 GeV; results from 2.76 TeV data for the 2-3 GeV interval were 
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obtained using wider intervals to improve statistics. Results are shown in the top panels of Fig. 4; 
the left and right panels show 2.76 and 13 TeV data, respectively. A significant dependence is seen. 
Separately, the same analysis was applied requiring both and Pj to fall within the above intervals. If 
factorization holds, the V 2 values calculated using: 

VliPTl) = V2,2{P1\,P12)I V^2,2(FT2>FT2)> (6) 

where pji and pji indicate which of the three intervals, 0.5-5, 0.5-1, and 2-3 GeV, p^ and p!j. are 
required to lie within, should be independent of pT 2 - The V 2 values obtained using Eq. 6 are shown in the 
middle panels of Fig. 4. For both collision energies, the three sets of V 2 values agree within uncertainties, 
indicating that V 22 factorizes. 

This analysis is sensitive to potential -dependent changes in the shape of the peripheral reference. For 
example, the pythia 8 sample shows a modest -dependent change in the width of the dijet peak for 
small Also, the V 2 p could vary with A^“ over the 0<AJ^‘^<20 range. To test the sensitivity of the 
results presented here to such shape changes, the analysis was repeated using 0-5, 0-10, and 10-20 A^®‘^ 
intervals to form TP“‘p’’(A 0). The largest resulting change in V 22 was taken as a systematic uncertainty. 
The relative uncertainty varies from 6% at A^“=30 to 2% for AJ^^^^hO in the 13 TeV data, and is less 
than <6% for all A^“ for the 2.76 TeV data. When using the 0-5 interval for TP'^"p*’(A 0), t; 2,2 values 
consistent with those shown in Fig. 4 are measured in A^^^^ intervals 5-10, 10-15 and 15-20. 

Potential systematic uncertainties on V 22 due to a residual l^cf) dependence of the two-particle acceptance 
that does not cancel in the SIB ratio are evaluated following Ref. [47] and are found to be less than 1%. 
The effect of the uncertainty on the tracking efficiency on 1 ) 2,2 is determined to be less than 1%. A separate 
systematic on V 22 due to the 0 and pj resolution of the charged-particle measurement is estimated to be 
2% (6%) for pt>0.5 GeV (pt<0.5 GeV). Events with unresolved multiple vertices decrease the measured 
1 ) 2,2 by increasing the combinatoric pedestal in T(A0) without increasing the modulation. The resulting 
systematic on V 22 increases with and is estimated to be less than 0.25% and 5% for the 13 and 
2.76 TeV data, respectively. The combined systematic uncertainties on V 2,2 and on V 2 are shown by the 
shaded boxes in Fig. 4. The total V 22 systematic uncertainty for 0.5<p^’^<5 GeV varies between ~5% at 
low A'^^^ to ~3% at high A^“ in the 13 TeV data, while in the 2.76 TeV data the uncertainty is 8% for all 
A^“. The systematic uncertainty on V 2 is approximately half that for ^ 2 , 2 - 

As shown in Fig. 4, the measured V 2 are independent of and are consistent between the two collision 
energies within uncertainties. The px dependence of V 2 for the 50-60 A^^‘^ interval, shown in the bottom 
left panel of Fig. 4, is similar for both collision energies to that previously measured in p-i-Pb and Pb-tPb 
collisions. It increases with px at low px, reaches a maximum between 2 and 3 GeV, and then decreases 
at higher px. The bottom right panel of Fig. 4 shows the px dependence of V 2 for different A^“ intervals; 
no significant dependence is observed. 

In summary, ATFAS has measured the multiplicity and px dependence of two-charged-particle corre¬ 
lations in V^=13 and 2.76 TeV pp collisions at the FHC. The correlation functions at both energies 
show a ridge whose strength increases with multiplicity. A new template fitting procedure shows that the 
per-trigger-particle yields for |A7/|>2 are described well by a superposition of the yields measured in a 
low-multiplicity interval and a constant modulated by cos (2A0). Thus, as observed in p-i-Pb collisions 
[4], the pp data presented here are compatible with both a “near-side” ridge centered at Ad> = 0 and an 
“away-side” ridge centered at Ad> = tt that both result from a sinusoidal component of the two-particle 
correlation. The extracted Fourier coefficients, 1 ) 2 , 2 , exhibit factorization, which is characteristic of a 
global modulation of the per-event single-particle distributions also seen in p-i-Pb and Pb-i-Pb collisions. 
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Figure 4; Measured 112,2 (top) and V 2 (middle) values versus for different intervals for 2.76 (left) and 13 TeV 
(right) data. Results are averaged over bins of width 10 spanning the range 20<AtJ^‘^<100 and 20<AtJ“<130 
for 2.76 and 13 TeV data, respectively, except for the 2<plj,<3 GeV results for the 2.76 TeV data which are averaged 
over bins of width 20. Measured V 2 values versus Pj (bottom) spanning the range 0.3</7iJ,<5.0 GeV for 13 and 
2.76 TeV data for the 50<V‘“<60 interval (left) and for three intervals in the 13 TeV data (right). Results 
are averaged over the intervals indicated by horizontal error bars. On all points, the vertical error bars indicate 
statistical uncertainties. The shaded bands indicate systematic uncertainties. For clarity, they are only shown for 
the 0.5<p)j,<5 GeV case in the middle, for 2.76 TeV data in the lower left, and for the 40<VJ“<50 case in the lower 
right panels. 

















The amplitudes, V 2 , of the single-particle modulation, are -independent and agree between 2.76 and 
13 TeV within uncertainties. They increase with px for px ^ 3 GeV and then decrease at higher px, 
following a trend similar to that observed in p-i-Pb and Pb-i-Pb collisions. These results suggest that the 
ridges in pp and p-i-Pb collisions may arise from a similar physical mechanism which does not have a 
strong yfs dependence. 
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